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SUMMARY 


i • t, ^ hree b °° ks on supercritical wing sections have been pub- 
■ rec ® nbly tha * document and list codes for the design and 
analysis of transonic airfoils. The codes have had a significant 
impact on the development of supercritical wing technologv . This 
li dev ° ted to several new contributions to the theory on 
which the cedes are based. It has become possible to prescribe 
P ressure distribution within a reasonable tolerance even over 
the supersonic portion of a shockless airfoil. The purely sub- 
sonic probiem of modeling the flow near the trailing edge has 
been handled in a way that eliminates any appreciable loss of 
lift in practice. Boundary layer effects are taken into account 
in an empirically satisfactory fashion. The work has been ex- 
en e to the case of cascades of airfoils appropriate for the 
design of compressor and turbine blades. These methods of com- 
putational fluid dynamics have produced a family of airfoils that 
deliver outstanding performance over a wide range of conditions. 


INTRODUCTION 


In the last few years three books have appeared that list 
computer codes for the design and analysis of transonic airfoils 
’ 4 . •' ?• T he design code relies on the method of complex 

characteristics in the hodograph plane to construct shockless 
airfoixs. The analysis code uses artifical viscosity to calcu- 
late flows with weak shock waves at off-design conditions. Com- 
parisons with experiments show that an excellent simulation of 
two-dimensional wind tunnel tests is obtained. The codes have 
been widely adopted by the aircraft industry as a tool for the 
development of supercritical wing technology. 


SYMBOLS 

f complex analytic function 

h real function 

*Work supportedljy NASA under Grants NGR-33-016-167 and NGR- 
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i square root of -1 

k real constant 

M Mach number 

q speed 

u function of x and y 

even function of y 
a 2 even function of y 

x canonical coordinate 

y canonical coordinate 

K characteristic coordinate 

h characteristic coordinate 

0 flow angle 

p density 

<J> velocity potential 

stream function 


ANALYSIS AND DESIGN 


For analysis, differencing of the partial differential 
equations of gas dynamics that do not adhere to strict con- 
servation form turns out to give the best representation of 
boundary layer-shock wave interaction over a broad range of 
conditions (ref . 2) . An improved formula for the wave drag 

compensates for errors in the conservation of mass across shocks 
that are of the third order in the shock strength. Modern tech- 
niques of conformal mapping and fast Fourier transform have led 
to an upgraded analysis code that has unusual speed and accuracy 

Corrections for the displacement thickness of the turbulent 
boundary layer have been made in both the design and the analysis 
codes. Even more important is an adequate model of the flow near 
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the trailing edge of the airfoil. It suffices to represent the 
wake by a pair of parallel streamlines across which the pressure 
balances. Large favorable pressure gradients can be tolerated on 
the lower surface to provide for heavy aft loading. However, a 
Stratford distribution should be used at the rear of the upper 
surface to limit the adverse pressure gradient there so as to 
avoid separation. Airfoils designed with this in mind perform 
well over a wide range of conditions. 

Shockless airfoils serve as an acceptable mathematical 
model for the design of supercritical wing sections. Drag creep 
can be reduced by restricting the size of the supersonic zone in 
the shockless flow. The design method has been extended to in- 
clude cascades of airfoils such as occur in compressors and tur- 
bines . Present codes can handle gap-to-cnord ratios down to 
unity (see figures 1 and 2) . The concept of a supercritical com- 
pressor blade has been tested successfully by Harry Stephens of 
Pratt and Whitney Aircraft in a cascade wind tunnel of the DFVLR 
in Germany (see figure 3) . 

The latest version of the design code enables one to assign 
the pressure distribution with a certain tolerance and still ob- 
tain shockless flow when it exists (ref. 3) . This is achieved 
by formulating a new boundary value problem in the unit circle 
of the complex plane of one of the characteristic coordinates for 
the gas dynamics equations. It is worthwhile to review briefly 
the theory underlying the new code, which seems to be quite suc- 
cessful in practice. 


In terms of characteristic coordinates ? and n , the partial 
differential equations for the velocity potential <J> and stream 
function ip of plane compressible flow can be expressed in the form 


= i 


/l-M 2 
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where M is the local Mach number and p is the density. Through 
analytic continuation the equations remain valid in the complex 
domain. The speed q and angle 0 of the flow are related to ? 
and n by the formulas 


log f (?) = | /l-M 2 I 3 - - i0 
log f (n) =| /lH? ^ + i6 , 

where f is an analytic function mapping the flow onto a region 
that can be taken as the unit circle |?| < 1. Analytic contin- 
uation around the sonic line can be performed along paths on 


which 1-M 2 does not vanish. 

The problem of finding an ^rSe^rob-* 

prescribed function °* ‘ h ® a 3hich the speed is assigned as a 
f^tL^"^)^ the velocity Potential 

iL% 1 aS i ? h rf^u^tirfrrrL a undar y conditions of the 
form 

Rs{ log f(U> = h(q) 

Re{ ij/ ( £ , 5 ) _ ik <f» ( 5 • C) ) = ® • 

where h is known in terms ?ri£l«d\te«tiveiy 

This nonlinear boundary valu P calcu i a ted and then com- 

by first guessing i so that f can be cai.ui For an ap - 

propriate a chot ce°of h h^and^k^the^iterations^converge^even^in^t^e^^ 

S’sSasSreri ansrrs.'K - “ 


u xx + u yy + 


u =0 

y 


obtained by bringing the Tricomi Ration /“ifctr c^with 

S^SSSTvSSrSfWoS ^boundary. Sant introduces a 
decomposition 

u(x,y) = y 2/ \(*.y 2 l + “2 ,x ' y2) 

1'is *”£ is'^led'to'boundary SSuTSSl^^l’p^ 
unit 1 semicircle that can be solved in closed form. 

Numerous examples of shockless airfoils 
by the method of complex characteristics g h two super- 

scribed qW). These include =^“® r a “ r£oils suitable for 
sonic zones (figs. 4 and ’• airfoils with large turning angles 
(fig? r 4 ). (f The procedure U is largely automatic, and it is rela- 
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tively easy to implement new ideas about design that depend on 
the pressure distribution. However, further conformal trans- 
formation of the characteristic coordinates will be required if 
9 a P“ to “Chord ratios significantly lower than unity are desired. 

It would be of interest to generalize the method of designing 
shockless airfoils based on giving a pressure distribution to the 
case of three-dimensional transonic flow past a swept wing. This 
might be achieved by modifying the analysis code appropriately 
and introducing an artifical viscosity that smears shocks. To 
start with it would be helpful to have a fast STAR version of the 
swept wing code published recently by Jameson and Caughey (refs. 

2 and 5) . Work on these proposals is in progress at the Courant 
Mathematics and Computing Laboratory. 
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M2-.762 DEL THr 99.26 


+ OUTPUT 


G/C=l .10 


Figure 2.- Sanz turbine blade. 



PRATT AND WHITNEY COMPRESSOR AIRFOIL Rrl.l MILLION 

— THEORY MU. 780 M2=.480 DEL TH=25.0 

A EXPERIMENT Ml=.775 M2r.544 DEL TH=25.5 L0SS=.0196 

Figure 3.- Supercritical compressor test data. 
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